Genetic underpinnings of arthropod community distributions in Populus trichocarpa.

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Ectomycorrhizal symbiosis prepares its host locally and systemically for abiotic cue signaling.

Tree growth and survival are dependent on their ability to perceive signals, integrate them, and trigger timely and fitted molecular and growth responses. While ectomycorrhizal symbiosis is a predominant tree-microbe interaction in forest ecosystems, little is known about how and to what extent it helps trees cope with environmental changes. We hypothesized that the presence of Laccaria bicolor influences abiotic cue perception by Populus trichocarpa and the ensuing signaling cascade. We submitted ectomycorrhizal or non-ectomycorrhizal P. trichocarpa cuttings to short-term cessation of watering or ozone fumigation to focus on signaling networks before the onset of any physiological damage. Poplar gene expression, metabolite levels, and hormone levels were measured in several organs (roots, leaves, mycorrhizas) and integrated into networks. We discriminated the signal responses modified or maintained by ectomycorrhization. Ectomycorrhizas buffered hormonal changes in response to short-term environmental variations systemically prepared the root system for further fungal colonization and alleviated part of the root abscisic acid (ABA) signaling. The presence of ectomycorrhizas in the roots also modified the leaf multi-omics landscape and ozone responses, most likely through rewiring of the molecular drivers of photosynthesis and the calcium signaling pathway. In conclusion, P. trichocarpa-L. bicolor symbiosis results in a systemic remodeling of the host's signaling networks in response to abiotic changes. In addition, ectomycorrhizal, hormonal, metabolic, and transcriptomic blueprints are maintained in response to abiotic cues, suggesting that ectomycorrhizas are less responsive than non-mycorrhizal roots to abiotic challenges.

Chromosome-level genome assembly of Quercus variabilis provides insights into the molecular mechanism of cork thickness.

Quercus variabilis is a deciduous woody species with high ecological and economic value, and is a major source of cork in East Asia. Cork from thick softwood sheets have higher commercial value than those from thin sheets. It is extremely difficult to genetically improve Q. variabilis to produce high quality softwood due to the lack of genomic information. Here, we present a high-quality chromosomal genome assembly for Q. variabilis with length of 791,89 Mb and 54,606 predicted genes. Comparative analysis of protein sequences of Q. variabilis with 11 other species revealed that specific and expanded gene families were significantly enriched in the "fatty acid biosynthesis" pathway in Q. variabilis, which may contribute to the formation of its unique cork. Based on weighted correlation network analysis of time-course (i.e., five important developmental ages) gene expression data in thick-cork versus thin-cork genotypes of Q. variabilis, we identified one co-expression gene module associated with the thick-cork trait. Within this co-expression gene module, 10 hub genes were associated with suberin biosynthesis. Furthermore, we identified a total of 198 suberin biosynthesis-related new candidate genes that were up-regulated in trees with a thick cork layer relative to those with a thin cork layer. Also, we found that some genes related to cell expansion and cell division were highly expressed in trees with a thick cork layer. Collectively, our results revealed that two metabolic pathways (i.e., suberin biosynthesis, fatty acid biosynthesis), along with other genes involved in cell expansion, cell division, and transcriptional regulation, were associated with the thick-cork trait in Q. variabilis, providing insights into the molecular basis of cork development and knowledge for informing genetic improvement of cork thickness in Q. variabilis and closely related species.

Validation of a metabolite-GWAS network for Populus trichocarpa family 1 UDP-glycosyltransferases.

Metabolite genome-wide association studies (mGWASs) are increasingly used to discover the genetic basis of target phenotypes in plants such as Populus trichocarpa, a biofuel feedstock and model woody plant species. Despite their growing importance in plant genetics and metabolomics, few mGWASs are experimentally validated. Here, we present a functional genomics workflow for validating mGWAS-predicted enzyme-substrate relationships. We focus on uridine diphosphate-glycosyltransferases (UGTs), a large family of enzymes that catalyze sugar transfer to a variety of plant secondary metabolites involved in defense, signaling, and lignification. Glycosylation influences physiological roles, localization within cells and tissues, and metabolic fates of these metabolites. UGTs have substantially expanded in P. trichocarpa, presenting a challenge for large-scale characterization. Using a high-throughput assay, we produced substrate acceptance profiles for 40 previously uncharacterized candidate enzymes. Assays confirmed 10 of 13 leaf mGWAS associations, and a focused metabolite screen demonstrated varying levels of substrate specificity among UGTs. A substrate binding model case study of UGT-23 rationalized observed enzyme activities and mGWAS associations, including glycosylation of trichocarpinene to produce trichocarpin, a major higher-order salicylate in P. trichocarpa. We identified UGTs putatively involved in lignan, flavonoid, salicylate, and phytohormone metabolism, with potential implications for cell wall biosynthesis, nitrogen uptake, and biotic and abiotic stress response that determine sustainable biomass crop production. Our results provide new support for in silico analyses and evidence-based guidance for in vivo functional characterization.

Expression quantitative trait loci mapping identified PtrXB38 as a key hub gene in adventitious root development in Populus.

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

Cell wall response of field grown Populus to Septoria infection.

Due to its ability to spread quickly and result in tree mortality, Sphaerulina musiva (Septoria) is one of the most severe diseases impacting Populus. Previous studies have identified that Septoria infection induces differential expression of phenylpropanoid biosynthesis genes. However, more extensive characterization of changes to lignin in response to Septoria infection is lacking. To study the changes of lignin due to Septoria infection, four field grown, naturally variant Populus trichocarpa exhibiting visible signs of Septoria infection were sampled at health, infected, and reaction zone regions for cell wall characterization. Fourier transform infrared spectroscopy (FTIR), nuclear magnetic resonance (NMR), and acid hydrolysis were applied to identify changes to the cell wall, and especially lignin. FTIR and subsequent principal component analysis revealed that infected and reaction zone regions were similar and could be distinguished from the non-infected (healthy) region. NMR results indicated the general trend that infected region had a higher syringyl:guaiacyl ratio and lower p-hydroxybenzoate content than the healthy regions from the same genotype. Finally, Klason lignin content in the infected and/or reaction zone regions was shown to be higher than healthy region, which is consistent with previous observations of periderm development and metabolite profiling. These results provide insights on the response of Populus wood characteristics to Septoria infection, especially between healthy and infected region within the same genotype.

Novel candidate genes for lignin structure identified through genome-wide association study of naturally varying Populus trichocarpa.

Populus is a promising lignocellulosic feedstock for biofuels and bioproducts. However, the cell wall biopolymer lignin is a major barrier in conversion of biomass to biofuels. To investigate the variability and underlying genetic basis of the complex structure of lignin, a population of 409 three-year-old, naturally varying Populus trichocarpa genotypes were characterized by heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR). A subsequent genome-wide association study (GWAS) was conducted using approximately 8.3 million single nucleotide polymorphisms (SNPs), which identified 756 genes that were significantly associated (-log10(p-value)>6) with at least one lignin phenotype. Several promising candidate genes were identified, many of which have not previously been reported to be associated with lignin or cell wall biosynthesis. These results provide a resource for gaining insights into the molecular mechanisms of lignin biosynthesis and new targets for future genetic improvement in poplar.

Diel dynamics of multi-omics in elkhorn fern provide new insights into weak CAM photosynthesis.

Crassulacean acid metabolism (CAM) has high water-use efficiency (WUE) and is widely recognized to have evolved from C3 photosynthesis. Different plant lineages have convergently evolved CAM, but the molecular mechanism that underlies C3-to-CAM evolution remains to be clarified. Platycerium bifurcatum (elkhorn fern) provides an opportunity to study the molecular changes underlying the transition from C3 to CAM photosynthesis because both modes of photosynthesis occur in this species, with sporotrophophyll leaves (SLs) and cover leaves (CLs) performing C3 and weak CAM photosynthesis, respectively. Here, we report that the physiological and biochemical attributes of CAM in weak CAM-performing CLs differed from those in strong CAM species. We investigated the diel dynamics of the metabolome, proteome, and transcriptome in these dimorphic leaves within the same genetic background and under identical environmental conditions. We found that multi-omic diel dynamics in P. bifurcatum exhibit both tissue and diel effects. Our analysis revealed temporal rewiring of biochemistry relevant to the energy-producing pathway (TCA cycle), CAM pathway, and stomatal movement in CLs compared with SLs. We also confirmed that PHOSPHOENOLPYRUVATE CARBOXYLASE KINASE (PPCK) exhibits convergence in gene expression among highly divergent CAM lineages. Gene regulatory network analysis identified candidate transcription factors regulating the CAM pathway and stomatal movement. Taken together, our results provide new insights into weak CAM photosynthesis and new avenues for CAM bioengineering.

Rapid screening of secondary aromatic metabolites in Populus trichocarpa leaves.

BACKGROUND: High-throughput metabolomics analytical methodology is needed for population-scale studies of bioenergy-relevant feedstocks such as poplar (Populus sp.). Here, the authors report the relative abundance of extractable aromatic metabolites in Populus trichocarpa leaves rapidly estimated using pyrolysis-molecular beam mass spectrometry (py-MBMS). Poplar leaves were analyzed in conjunction with and validated by GC/MS analysis of extracts to determine key spectral features used to build PLS models to predict the relative composition of extractable aromatic metabolites in whole poplar leaves. RESULTS: The Pearson correlation coefficient for the relative abundance of extractable aromatic metabolites based on ranking between GC/MS analysis and py-MBMS analysis of the Boardman leaf set was 0.86 with R2 = 0.76 using a simplified prediction approach from select ions in MBMS spectra. Metabolites most influential to py-MBMS spectral features in the Clatskanie set included the following compounds: catechol, salicortin, salicyloyl-coumaroyl-glucoside conjugates, α-salicyloylsalicin, tremulacin, as well as other salicylates, trichocarpin, salicylic acid, and various tremuloidin conjugates. Ions in py-MBMS spectra with the highest correlation to the abundance of extractable aromatic metabolites as determined by GC/MS analysis of extracts, included m/z 68, 71, 77, 91, 94, 105, 107, 108, and 122, and were used to develop the simplified prediction approach without PLS models or a priori measurements. CONCLUSIONS: The simplified py-MBMS method is capable of rapidly screening leaf tissue for relative abundance of extractable aromatic secondary metabolites to enable prioritization of samples in large populations requiring comprehensive metabolomics that will ultimately inform plant systems biology models and advance the development of optimized biomass feedstocks for renewable fuels and chemicals.

Mechanism for Utilization of the Populus-Derived Metabolite Salicin by a Pseudomonas-Rahnella Co-Culture.

Pseudomonas fluorescens GM16 associates with Populus, a model plant in biofuel production. Populus releases abundant phenolic glycosides such as salicin, but P. fluorescens GM16 cannot utilize salicin, whereas Pseudomonas strains are known to utilize compounds similar to the aglycone moiety of salicin-salicyl alcohol. We propose that the association of Pseudomonas to Populus is mediated by another organism (such as Rahnella aquatilis OV744) that degrades the glucosyl group of salicin. In this study, we demonstrate that in the Rahnella-Pseudomonas salicin co-culture model, Rahnella grows by degrading salicin to glucose 6-phosphate and salicyl alcohol which is secreted out and is subsequently utilized by P. fluorescens GM16 for its growth. Using various quantitative approaches, we elucidate the individual pathways for salicin and salicyl alcohol metabolism present in Rahnella and Pseudomonas, respectively. Furthermore, we were able to establish that the salicyl alcohol cross-feeding interaction between the two strains on salicin medium is carried out through the combination of their respective individual pathways. The research presents one of the potential advantages of salicyl alcohol release by strains such as Rahnella, and how phenolic glycosides could be involved in attracting multiple types of bacteria into the Populus microbiome.

Exploring the role of plant lysin motif receptor-like kinases in regulating plant-microbe interactions in the bioenergy crop Populus.

For plants, distinguishing between mutualistic and pathogenic microbes is a matter of survival. All microbes contain microbe-associated molecular patterns (MAMPs) that are perceived by plant pattern recognition receptors (PRRs). Lysin motif receptor-like kinases (LysM-RLKs) are PRRs attuned for binding and triggering a response to specific MAMPs, including chitin oligomers (COs) in fungi, lipo-chitooligosaccharides (LCOs), which are produced by mycorrhizal fungi and nitrogen-fixing rhizobial bacteria, and peptidoglycan in bacteria. The identification and characterization of LysM-RLKs in candidate bioenergy crops including Populus are limited compared to other model plant species, thus inhibiting our ability to both understand and engineer microbe-mediated gains in plant productivity. As such, we performed a sequence analysis of LysM-RLKs in the Populus genome and predicted their function based on phylogenetic analysis with known LysM-RLKs. Then, using predictive models, molecular dynamics simulations, and comparative structural analysis with previously characterized CO and LCO plant receptors, we identified probable ligand-binding sites in Populus LysM-RLKs. Using several machine learning models, we predicted remarkably consistent binding affinity rankings of Populus proteins to CO. In addition, we used a modified Random Walk with Restart network-topology based approach to identify a subset of Populus LysM-RLKs that are functionally related and propose a corresponding signal transduction cascade. Our findings provide the first look into the role of LysM-RLKs in Populus-microbe interactions and establish a crucial jumping-off point for future research efforts to understand specificity and redundancy in microbial perception mechanisms.

Mechanistic Insights into Cell-Free Gene Expression through an Integrated -Omics Analysis of Extract Processing Methods.

Cell-free systems derived from crude cell extracts have developed into tools for gene expression, with applications in prototyping, biosensing, and protein production. Key to the development of these systems is optimization of cell extract preparation methods. However, the applied nature of these optimizations often limits investigation into the complex nature of the extracts themselves, which contain thousands of proteins and reaction networks with hundreds of metabolites. Here, we sought to uncover the black box of proteins and metabolites in Escherichia coli cell-free reactions based on different extract preparation methods. We assess changes in transcription and translation activity from σ70 promoters in extracts prepared with acetate or glutamate buffer and the common post-lysis processing steps of a runoff incubation and dialysis. We then utilize proteomic and metabolomic analyses to uncover potential mechanisms behind these changes in gene expression, highlighting the impact of cold shock-like proteins and the role of buffer composition.

Multi-omic characterization of bifunctional peroxidase 4-coumarate 3-hydroxylase knockdown in Brachypodium distachyon provides insights into lignin modification-associated pleiotropic effects.

A bifunctional peroxidase enzyme, 4-coumarate 3-hydroxylase (C3H/APX), provides a parallel route to the shikimate shunt pathway for the conversion of 4-coumarate to caffeate in the early steps of lignin biosynthesis. Knockdown of C3H/APX (C3H/APX-KD) expression has been shown to reduce the lignin content in Brachypodium distachyon. However, like many other lignin-modified plants, C3H/APX-KDs show unpredictable pleiotropic phenotypes, including stunted growth, delayed senescence, and reduced seed yield. A system-wide level understanding of altered biological processes in lignin-modified plants can help pinpoint the lignin-modification associated growth defects to benefit future studies aiming to negate the yield penalty. Here, a multi-omic approach was used to characterize molecular changes resulting from C3H/APX-KD associated lignin modification and negative growth phenotype in Brachypodium distachyon. Our findings demonstrate that C3H/APX knockdown in Brachypodium stems substantially alters the abundance of enzymes implicated in the phenylpropanoid biosynthetic pathway and disrupt cellular redox homeostasis. Moreover, it elicits plant defense responses associated with intracellular kinases and phytohormone-based signaling to facilitate growth-defense trade-offs. A deeper understanding along with potential targets to mitigate the pleiotropic phenotypes identified in this study could aid to increase the economic feasibility of lignocellulosic biofuel production.

Comparative genomics analysis of drought response between obligate CAM and C3 photosynthesis plants.

Crassulacean acid metabolism (CAM) plants exhibit elevated drought and heat tolerance compared to C3 and C4 plants through an inverted pattern of day/night stomatal closure and opening for CO2 assimilation. However, the molecular responses to water-deficit conditions remain unclear in obligate CAM species. In this study, we presented genome-wide transcription sequencing analysis using leaf samples of an obligate CAM species Kalanchoë fedtschenkoi under moderate and severe drought treatments at two-time points of dawn (2-h before the start of light period) and dusk (2-h before the dark period). Differentially expressed genes were identified in response to environmental drought stress and a whole genome wide co-expression network was created as well. We found that the expression of CAM-related genes was not regulated by drought stimuli in K. fedtschenkoi. Our comparative analysis revealed that CAM species (K. fedtschenkoi) and C3 species (Arabidopsis thaliana, Populus deltoides 'WV94') share some common transcriptional changes in genes involved in multiple biological processes in response to drought stress, including ABA signaling and biosynthesis of secondary metabolites.

Relationships between Sphaerulina musiva Infection and the Populus Microbiome and Metabolome.

Pathogenic fungal infections in plants may, in some cases, lead to downstream systematic impacts on the plant metabolome and microbiome that may either alleviate or exacerbate the effects of the fungal pathogen. While Sphaerulina musiva is a well-characterized fungal pathogen which infects Populus tree species, an important wood fiber and biofuel feedstock, little is known about its systematic effects on the metabolome and microbiome of Populus. Here, we investigated the metabolome of Populus trichocarpa and Populus deltoides leaves and roots and the microbiome of the leaf and root endospheres, phylloplane, and rhizosphere to understand the systematic impacts of S. musiva abundance and infection on Populus species in a common garden field setting. We found that S. musiva is indeed present in both P. deltoides and P. trichocarpa, but S. musiva abundance was not statistically related to stem canker onset. We also found that the leaf and root metabolomes significantly differ between the two Populus species and that certain leaf metabolites, particularly the phenolic glycosides salirepin and salireposide, are diminished in canker-infected P. trichocarpa trees compared to their uninfected counterparts. Furthermore, we found significant associations between the metabolome, S. musiva abundance, and microbiome composition and α-diversity, particularly in P. trichocarpa leaves. Our results show that S. musiva colonizes both resistant and susceptible hosts and that the effects of S. musiva on susceptible trees are not confined to the site of canker infection. IMPORTANCE Poplar (Populus spp.) trees are ecologically and economically important trees throughout North America. However, many western North American poplar plantations are at risk due to the introduction of the nonnative fungal pathogen Sphaerulina musiva, which causes leaf spot and cankers, limiting their production. To better understand the interactions among the pathogen S. musiva, the poplar metabolome, and the poplar microbiome, we collected leaf, root, and rhizosphere samples from poplar trees consisting of 10 genotypes and two species with differential resistance to S. musiva in a common garden experiment. Here, we outline the nuanced relationships between the poplar metabolome, microbiome, and S. musiva, showing that S. musiva may affect poplar trees in tissues distal to the site of infection (i.e., stem). Our research contributes to improving the fundamental understanding of S. musiva and Populus sp. ecology and the utility of a holobiont approach in understanding plant disease.

Reporter genes confer new-to-nature ornamental traits in plants.

This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).

Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon.

Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA.

As the focus for CRISPR/Cas-edited plants moves from proof-of-concept to real-world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one guide RNA (gRNA) per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in hybrid poplar (Populus tremula × P. alba INRA 717-1B4). Unexpected duplications of MYB186 and MYB138 resulted in eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and polymerase chain reaction analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole-leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the nonglandular trichomes of poplar.

Developmental changes in lignin composition are driven by both monolignol supply and laccase specificity.

The factors controlling lignin composition remain unclear. Catechyl (C)-lignin is a homopolymer of caffeyl alcohol with unique properties as a biomaterial and precursor of industrial chemicals. The lignin synthesized in the seed coat of Cleome hassleriana switches from guaiacyl (G)- to C-lignin at around 12 to 14 days after pollination (DAP), associated with a rerouting of the monolignol pathway. Lack of synthesis of caffeyl alcohol limits C-lignin formation before around 12 DAP, but coniferyl alcohol is still synthesized and highly accumulated after 14 DAP. We propose a model in which, during C-lignin biosynthesis, caffeyl alcohol noncompetitively inhibits oxidation of coniferyl alcohol by cell wall laccases, a process that might limit movement of coniferyl alcohol to the apoplast. Developmental changes in both substrate availability and laccase specificity together account for the metabolic fates of G- and C-monolignols in the Cleome seed coat.

Diversity and conservation of plant small secreted proteins associated with arbuscular mycorrhizal symbiosis.

Arbuscular mycorrhizal symbiosis (AMS) is widespread mutualistic association between plants and fungi, which plays an essential role in nutrient exchange, enhancement in plant stress resistance, development of host, and ecosystem sustainability. Previous studies have shown that plant small secreted proteins (SSPs) are involved in beneficial symbiotic interactions. However, the role of SSPs in the evolution of AMS has not been well studied yet. In this study, we performed computational analysis of SSPs in 60 plant species and identified three AMS-specific ortholog groups containing SSPs only from at least 30% of the AMS species in this study and three AMS-preferential ortholog groups containing SSPs from both AMS and non-AMS species, with AMS species containing significantly more SSPs than non-AMS species. We found that independent lineages of monocot and eudicot plants contained genes in the AMS-specific ortholog groups and had significant expansion in the AMS-preferential ortholog groups. Also, two AMS-preferential ortholog groups showed convergent changes, between monocot and eudicot species, in gene expression in response to arbuscular mycorrhizal fungus Rhizophagus irregularis. Furthermore, conserved cis-elements were identified in the promoter regions of the genes showing convergent gene expression. We found that the SSPs, and their closely related homologs, in each of three AMS-preferential ortholog groups, had some local variations in the protein structural alignment. We also identified genes co-expressed with the Populus trichocarpa SSP genes in the AMS-preferential ortholog groups. This first plant kingdom-wide analysis on SSP provides insights on plant-AMS convergent evolution with specific SSP gene expression and local diversification of protein structures.